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INTRODUCTION 

Historically  the  problem  of  propagation  of  optical  energy  through 
a scattering  senium  first  became  evident  in  astronomy.  More  recently 
interest  in  optical  communication  systems  employing  atmospheric  propa- 
gation channels  has  added  additional  impetus  to  the  search  for  analytic 
models  which  can  predict  system  performance. 

From  the  communications  system  viewpoint  one  may  ask:  How  does  a 

cloud  located  along  the  propagation  path  affect  the  communication  channel? 

A f miliar  observation  is  that  significant  amounts  of  sunlight  can  be 
observed  even  through  very  thick  clouds.  It  seems  reasonable  to  expect 
that  significant  amounts  of  optical  signal  could  be  collected  provided  the 
receiver's  f leld-of-view  (FOV)  is  increased  to  collect  much  of  the  scat- 
tered energy  in  addition  to  the  unscattered  component.  Such  a scheme 
introduces  multiple  paths  by  which  transmitted  photons  may  arrive  at  the 

receiver.  Thus  ^multipath  time  spreading,  of  the  optical  pulse  is  a 

5-  - T' 

* , j . 

consequence  which  must  be  accepted  when  receiver  FOV  accomodation  is  used! 

> 

Given  the  high  data  rates  expected  of  a viable  optical  communication 
system 'coupled  with  the  fact  that  the  information  is  transmitted  by  very 
short  (picosecond)  pulses,  mult ipa* h time  spreading  of  the  received 
pulses  can  produce  intersymbol  interference  increasing  the  error  rate.' 

Thus,  the  presence  of  a scattering  medium  in  the  propagation  path  may 
force  a redaction  of  the  data  rate  to  achieve  an  acceptable  error  rate. 

Although  numerous  theoretical  treatments  of  this  problem  can  be  found 

in  the  literature  [1,2,3]  very  few  experimental  results  have  been  published. 

* „ */~ 

”^Tbe  fall  scale  experiments  which  have  been  reported  typically  specify  the 
scattcrcr's  parameters  such  as  particle,  size  distribution,  concentration. 


2 

size  and  shape  ''with  little  precision.  Thus , one  objective  of  this  inves- 

? 

tlgation  was  to  conduct  a controlled  scattering  experiment  wherein  these 
parameters  were  known  quite  precisely.  The  uniform  latex  particles  in 
aqueous  suspension  which  were  used  to  simulate  water  droplets  in  the 
atmosphere  satisfied  this  requirement  well. 

A mode  locked  He:Ne  laser  provided  optical  pulses  typical  of  those 
encountered  in  optical  communications  systems.  JThe  time  delays  generated 
in  this  experiment  were  on  the  order  of  tens  of  picoseconds;  this  required 
an  indirect  measurement  scheme  which  compared  the  relative  phase  of  the 
incident  and  scattered  pulse  envelopes. 


SECTION  1 


THEORY 


A.  Mig  Solution 

Consider  a single  dielectric  sphere  whose  center  is  coincident  with 
the  origin  of  a spherical  coordinate  system  (Pig.  1-1).  Suppose  a 
uniform  plane  wave  propagating  in  the  +z-direction  is  incident  upon 
the  sphere  of  radius  a.  The  scattered  fields  are  given  by  the  solution 
usually  attributed  to  Mie  [4].  This  solution  is  exact  but  of  series  form. 

The  values  for  the  far  fields  (r  » X)  are  given  by: 


Eg  * ^ exp(-jkr]cos  <J>  S2(0) 


(1-1) 


(1-2) 


where  the  usual  exp  (jmt)  is  assumed  for  the  time  dependence  and 


S.(0)  = Z.  --/****  v {a  tt  (cost))  + b T (cose)} 
i n=JL  utn-rij  n n n n 


(1-3) 


S_(6)  = I.  -y—r  {b  * (cosO)  4-  a T (cosO)} 
2 n=l  n(n+l)  n n n n 


(1-4) 


(Y)  <{'  (8)  - m (y)  *'(3) 


n n 

««j^(Y)  ’}'n(S> 
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Where  a is  the  relative  refractive  index 


(1-8) 


^(z)  = sn(z>  (1-9) 

Cn<2)  = ^M+it  (z)  (1.10) 

•n  (c.os6)  * —r~a  P*(cos6)  (1-11) 

n s inti  n 

T (cos6)  « ~s  P1  (cosO)  (1.12) 

n do  n 


■* 

Equations  (1-9)  and  (1.10)  are  the  Riccat i -Bessel  functions  and  Px(cos0) 

n 

Is  an  associated  Legendre  polynoaial  in  equations  (1-11)  and  (1.12). 

Clearly  Lhe  computation  of  5^(8)  and  presents  a formidable 

calculation.  Fortunately  there  exist  tables  of  computed  values  for 
these  quantities.  If  only  the  intensity  of  the  field  is  of  interest  then 
j$^(Q)|^  and  1 (9 ) 1 ^ will  suffice. 

If  the  incident  field  is  linearly  polarized  then  the  scattered 
intensity  produced  by  a beam  polarized  perpendicular  to  the  observnt ion 
plane  is 

I --^j|S(0)|J  OU) 

An  r 


and  the  scattered  intensity  produced  by  a beam  polarized  parallel  to  the 
observation  plane  is 


1 

II 


,2  2 

An  r 


ls,(0)|2 


(1-1A) 
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The  computation  problem  then  reduces  to  evaluating  the  parameters  0 

and  a where 


2fla  c Particle  Diameter 
X Wavelength  in  Medium 


(1.15) 


*1  = Refractive  index  of  scatterer 
*2  Refractive  index  of  medium 


(1.16) 


2 2 

One  then  finds  the  tabulated  values  {s^(0)|  and  15,(0)!“  corresponding  to 
the  given  values  of  g and  m and  computes  the  intensity  of  the  scattered 
field  from  a single  sphere. 

Now  assume  an  ensemble  of  X identical,  randomly  distributed  spheres - 
tinder  this  assumption  the  scattered  fields  are  incoherent  at  any  arbitrary 
observation  point.  Thus  the  intensities  may  be  added  by  superposition  and 
the  intensity  from  the  ensemble  of  X scatterers  is  just  the  single  sphere 
intensity  multiplied  by  X. 

B.  Scattering  function  for  a single,  lossless  dielectric  sphere. 

Polystyrene  is  essentially  lossless  at  opvical  frequencies  so  that  the 
refractive  index  m^  = 1.5905  |5].  Ihe  refractive  index  m^  of  deionized 
water  is  = 1.332.  Therefore,  the  relative  refractive  index,  m of 
the  latex  particles  suspended  in  water  is 


m 


1.194 


(1.17) 


The  parameter 


B 


2na 

X 


3.181 


(1.18) 
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where  a is  the  particle  raJius  and  A is  the  optical  wavelength  in  the 


ieditua  having  refractive  index,  m^.  If  we  take 


6 = 3.2  = 3.181 


a - 1.2  = 1.194 


(1.19) 


2 2 

we  may  use  the  computed  values  for  |s^(0)|  and  | S2 (Q) | for  scattering 


angles  between  0°  and  180°  in  5°  steps  {6].  These  values  are  given  in 


Appendix  A. 


In  this  experiment  the  incident  beam  is  linearly  polarized  in  the  y 


direction.  The  observation  plane  is  the  x - z plane,  therefore,  the 


incident  bean  is  polarized  perpendicular  to  the  observation  plane. 


A polar  plot  of  1(9)  (Fig.  1-2)  illustrates  the  highly  forward 


character  of  the  scattered  intensity. 


If  a smooth  curve  is  used  to  connect  the  discrete  points,  an  analytic 


function  can  be  constructed  which  approximates  the  exact  Mie  series 


solution.  The  approximate  function,  F(9)  can  be  used  to  investigate 


multiple  scattering  analytically. 


A rectangular  plot  of  1(9)  (Fig.  1-3)  suggests  a Gaussian  function 


might  give  close  agreement  at  each  discrete  point,  let 


F(0)  = exp(-vt9*')  , -it  < 0 < 7i 


(1.20) 


Setting  F(0)  = 1(0)  at  the  half  amplitude  point  (0  = ±31°)  requires 


expl-n(.541)  ] = .5  , a = 2.368 


(1.21) 





F(0)  dO  = F(0)  dO 


w 

| exp(-2.368  O2)d0  » \ (^1*  = ^ 


a - , p<e)  - (l,!^,^-2-3680' 


A useful  quantity  which  can  be  calculated  now  is  the  average  cosine 
of  the  scattering  angle  6,  <cos  0> 


<cos  0>  = f cos  0 P(0)  d0  = 

-TT 


As  expected*  <cos  0>  indicates  a strong  forward  characteristic. 

C.  Multiple  versus  single  scattering 
1.  Analysis  of  single  scattering. 

Consider  the  detector  geometry  (Fig.  1-4).  The  following  definitions 
are  useful: 

<{>  - angular  displacement  of  detector 
20  - detector  f ield-of-view  (FOV) 

zq  - distance  from  incident  beam  entry  point  (z  = 0)  and  detector 
photocathode 

2(0.40  - intersection  of  incident  beam  with  limit  of  detector  FOV 

The  following  relations  obtain  from  Fig.  l-5a: 

a = z sin  0 (1.31] 

o 


m 


* - f - 0 - ♦ 


l£KUkMi>UI  fijtU 
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z(6,4>) 


z „ sin  0 
o 


cos  (J/  cos[^-(0  + <J>)] 


(1.33) 


z sin  0 

g cog  ^ + s^n  ^ CQS  g ZQ  C(®»^) 


(1.34) 


where 


C(M) 


sin  8 

sin  0 cos  $ + sin  <p  cos  8 


(1.35) 


For  0 < <j>  _<  12*  the  values  of  £(0,<J>)  and  z (0 »<f>)  (given  that  z = 31-75  cm 

o 

0 - 11#)  are: 


♦ (deg)  C(6»$)  z(8»<t>)  (cm) 


0 

1 

31.75 

1 

.918 

29.15 

2 

.848 

26.92 

3 

.789 

25.05 

4 

.737 

23.40 

5 

.692 

21.97 

6 

.653 

20.73 

7 

.617 

19.59 

8 

.586 

18.61 

9 

.558 

17.72 

10 

.532 

16.89 

11 

.509 

16.16 

12 

.488 

15.49 

Table  1-1  S(0,<{>)  and  z(0,<{>) 


Assume  a single  scattering  event  occurs  as  a given  photon  propagates 
through  the  chamber.  Let  2(6,1}))  be  the  coordinate  where  this  event  occurs. 
The  following  relations  obtain  from  Fig.  l-5b: 


s = z(0,<}>)  cos  4> 


(1.36) 


z - s z -•  z(0,A)  cos  <f> 
o o 


cos  6 


cos  0 


(1.37) 


z - z C(0,<}>)  cos  6 
o o 

cos  0 


(1.38) 


z (1  - £(0,<!>)  cos  (}») 

z*  - -° = 

cos  9 


(1.39) 


Note  that  the  change  in  path  length,  AL 


AL  = [z(0,$)  + z'l  - z 


(1.40) 


Equation  (1.40)  gives  the  limiting  value  of  AI.  for  single  scattering. 

The  time  delay.  At  for  a photon  traveling  the  differential  path  length  AT. 
is  given  by 


(1.41) 


where  c is  the  propagation  velocity  of  light  in  the  medium  having  refractive 


index  m^. 


Consider  the  special  case  of  propagating  optical  pulses  which  are 
periodic  with  period,  T.  Suppose  signals  are  derived  from  both  the  tin- 


— a 


scattered  main  beam  and  the  scattered  intensity  at  some  arbitrary  position. 

Of  course,  these  signals  retain  the  periodic  nature  of  the  optical  pulses 
and  a phase  difference  Aa  exists  between  the  two  periodic  signals.  Ue  note 
that 


Aa  At 
2tt  “ T 


(1.42) 


where  T = 13.14  nanoseconds.  Solving  Equations  (1.39)  through  (1.42)  we 
obtain  the  following  values: 


<|>(deg) 

z(0,4>)(cm) 

z’ (cm) 

0 

31.75 

0 

1 

29.16 

2.643 

2 

26.92 

4.937 

. 3 

25.06 

6.850 

4 

23.40 

8.564 

5 

21.98 

iQ.038 

6 

20.74 

11.332 

7 

19.60 

12.526 

8 

18.62 

13.560 

9 

17.72 

14.515 

10 

16.90 

15.389 

11 

16.16 

16.184 

12 

15.50 

16.899 

AL(cm) 

At(psec) 

An (deg) 

0 

• 0 

0 

.053 

2.35 

.064 

.107 

4.75 

.13 

.160 

7.10 

.20 

.214 

9.50 

.26 

.268 

11.88 

.33 

.322 

14.27 

.39 

.376 

16.68 

.46 

.430 

19.08 

.52 

.485 

21.49 

.59 

.539 

23.92 

.66 

.594 

26.35 

.72 

.649 

28.78 

.79 

Table  1-2 


iBw„ 
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Assume  now  that  the  single  scattering  event  occurs  at  z(0,<ft) . This 

2 

is  the  coordinate  of  the  intersection  of  the  incident  beam  with  the  re- 
stricted FOV,  9. 

The  following  relations  obtain  from  Fig-  l-5c; 


z(9,<fr) 


cos  ($1 


(1.43) 


2 - S 

O 


cos  2 


z(0,$) 

cos 


e 

cos  - 


* d - 
.«  _ o 


0 

COS  2 


AL  = + 2 * ) - 


(1.44) 


(1.45) 


(1.4*) 


Solving  equations  (1.45),  (1.46),  (1.41)  and  (1.42)  (given  that 

0 = 11*,  z ■ 31.75  cm)  we  obtain  the  values  shown  in  Table  1-3.  It  is 
o 

clear  from  the  results  summarized  in  Tables  1-2  and  1-3  that  a 
progressive  multipath  delay  results  at  viewing  positions  displaced 
from  the  main  beam.  However,  even  in  the  limiting  case  of  = 12°,  At 
is  less  than  30  picoseconds.  If  we  assume  a typical  pulse  width  of  300 
picoseconds,  this  delay  is  probably  negligible.  Indeed,  if  on  the 


average  only  a single  scattering  event  occurs,  then  much  of  the  signal 
energy  will  be  confined  to  the  main  beam  and  those  photons  which 


are  delayed  make  a negligible  contribution  to  the  total  received  pulse. 


Thus  multipath  time  spreading  from  single  scattering  is  negligible  In  so 


far  as  data  rate  limitation  is  concerned. 


z 


20 


[deg) 

*^<0 

z*  (cm) 

AL(cm) 

At(pscc) 

Aa(deg 

0 

15.89 

0 

0 

0 

0 

1 

14.58 

17.252 

.082 

3.62 

.10 

2 

13.46 

18.383 

.093 

4.12 

.11 

3 

12.53 

19.326 

.106 

4.71 

.13 

4 

11.70 

20.171 

.121 

5.38 

.15 

5 

10.99 

20.898 

.138 

6.12 

.17 

6 

10.37 

21.536 

.156 

6.91 

.19 

7 

9.80 

22.125 

.175 

7.75 

.21 

8 

9.31 

22.635 

.195 

8.64 

.24 

9 

8.86 

23.105 

.215 

9.55 

.26 

10 

8.45 

23.537 

.237 

10.50 

.29 

11 

8.05 

23.929 

.259 

11.47 

.31 

12 

7.75 

24.281 

.281 

12.46 

.34 

Table  1-3 


2.  Multiple  Scattering 

The  analysis  of  single  scattering  indicates  that  it  does  not  intro- 
duce significant  multipath  time  spreading.  Thus  multipath  scattering 
events  must  be  considered.  It  is  clear  that  a simple  geometrical  analysis 
analogous  to  that  used  for  the  single  scatter  case  does  not  exist. 

The  sultiple  scattering  problem  is  generally  analyzed  using  radiative 
transfer  theory  17]  or  by  performing  a Monte  Carlo  simulation  |S] . 

A typical  Monte  Carlo  photon  path  simulation  is  depicted  in  Fig.  1-6. 
The  incident  photon  travels  a distance  d^  to  the  i-tli  scatter  event.  A 
new  propagation  direction  is  selected  relative  to  the  i-tli  direction. 
Note  that  $ is  a function  of  the  scattering  angle.  0^  and  the  polar 


1 ~ ;'  1 


angle,  4>^-  The  process  is  repeated  until  the  photon  exits  the  simulated 
medium.  The  parameter  d^  is  a random  variable  with  probability  density 
P(d^)  related  to  the  mean  free  path  between  scatters.  0^  is  a random 
variable  with  probability  density  P(G^)  determined  by  scattering  function. 
4>£  is  a uniformly  distributed  random  variable.  Typically  the  total 
distance  traveled,  exit  angle  and  radius  of  the  exit  point  from  the  main 
beam  are  recorded  for  each  transmitted  photon.  After  a large  number 


of  these  events  have  been  collected,  the  data  is  reduced  to  obtain 
a representative  statistical  sample  of  all  paths.  An  analysis  of 
multipath  time  spreading  may  then  be  obtained  from  the  differential 
path  lengths. 
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SECTION  2 
EXPERIMENT 

A.  Equipment  and  configuration  description 

The  optical  source  is  a Spectra-Physics  Model  125  He:Ne  laser  operating 

© 

at  6328A  (Fig.  2-1).  The  mode  locked  output  is  a train  of  1.5  nanosecond 
(at  half  amplitude)  pulses  with  period  "f  13.14  nanoseconds  (76  MHz)  and 
average  power  of  25  milliwatts.  A beam  splitter  diverts  a small  amount 
of  the  beam  to  an  RCA  931A  photomultiplier  the  output  of  which  produces 
the  reference  RF  signal. 

The  main  beam  passes  through  a wedge  variable  attenuator  and  into 
the  scattering  chamber.  An  RCA  8644  Photomultiplier  detects  the  scattered 
intensity  at  any  desired  viewing  position.  Aperture  stops  control  the 
detector  f ield-of-view  and  a rotatable  polarizer  selects  the  desired 
orientation  of  polarization  of  scattered  intensity. 

Each  photomultiplier  output  is  followed  by  a narrow  pass-band 
filter  centered  at  the  fundamental  frequency,  thus  higher  order  frequency 
components  are  rejected. 

The  filtered  outputs  are  applied  to  the  sampling  probes  of  the 
Hewlett-Packard  8405A  Vector  Voltmeter.  The  RF  envelopes  are  frequency 
translated  to  an  intermediate  frequency  (IF)  and  the  IF  signals  are 
filtered  to  obtain  inputs  for  the  voltmeter  and  phasemeter. 

The  scattering  chamber  was  designed  so  that  the  geometric  path  length 
between  the  incident  beam  entry  point  and  the  viewing  position  at  the 
detector  was  independent  of  the  angled,  (Fig.  2-2).  This  was 
accomplished  by  making  the  output  window  of  the  chamber  the  surface  of 
a right  circular  cylinder  whose  center  is  located  at  the  incident 
beam  entry  point.  Tito  detector  is  mounted  on  a swing  arm  with  pivot 
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point  also  located  directly  above  the  incident  beam  entry  point. 

The  semi-infinite  extent  of  the  scattering  medium  (in  the  directions 
transverse  to  the  axis  of  propagation)  was  simulated  by  making  the  trans- 
verse dimensions  much  larger  than  the  diameter  of  the  incident  beam  and 
by  painting  the  interior  bottom,  side  and  top  surfaces  with  flat  black 
paint  to  minimize  reflections.  The  entrance  and  exit  windows  of  the 
chamber  were  fabricated  from  clear  plexiglass  acrylic  sheets.  The 
volume  of  the  chamber  is  approximately  3 liters. 

Consider  the  result  of  allowing  the  detector  to  scan  all  angles  4> 
such  that  the  received  signal  is  the  sum  over  all  <j>.  An  effective  field- 
of-view  can  be  defined  (Fig.  2-3). 


F°y  eff  « 26  + 2<j> 


The  effective  FOV  (given  that  0 = 11°,  4>  = 12°)  is  46°.  Such 
a FOV  is  unrealistically  large  for  a real  system  however  it  should 
be  noted  that  this  results  from  the  large  angles  needed  to  generate 
neasureable  multipath  delays  with  laboratory  scale  media. 

Initially  a highly  dilute  mixture  of  skim  milk  and  distilled 
water  was  used  as  a scattering  medium.  Casein  micelles  in  skim  milk 
act  as  highly  hydrated  spheres  of  mean  radius  800X.  Approximately  95% 

O 

of  these  are  in  the  range  40O-22OOA,  (9].  Thus  the  larger  micelles  have 

O 

strong  forward  scattering  characteristics  at  6328A.  Milk  has  the 
advantages  of  being  readily  available  and  relatively  inexpensive.  The 
disadvantages  are  short  shelf  life  and  the  uncertainty  as  to  exact 
distribution  of  particle  size  and  number. 
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Experiments  were  also  conducted  using  monodisperse  uniform  latex 
particles  manufactured  by  Dow  Diagnostics.  The  0.481  micron  diameter 
Polystyrene  particles  were  used  to  obtain  strong  forward  scattering.  These 
particles  are  spherical  while  suspended  in  deionized  water  and  have 
extremely  uniform  radii  (standard  deviation  18  A) . The  concentration 
(number  of  particles  per  milliliter)  was  calculated  from  technical  data 
supplied  by  the  manufacturer  and  verified  by  measuring  the  extinction 
coefficient  of  a dilute  suspension  of  the  particles.  (For  details  refer 
to  Appendix  B). 

The  concentrations  used  resulted  in  a scattering  medium  having  an 
optical  thickness  in  the  range  of  1-2.  (See  Appendix  C). 

B.  Measurement  procedure 

Initial  time  dispersion  measurements  were  attempted  by  observing  the 
scattered  pulse  in  real  time  using  the  Sylvania  502  crossed  field  photo- 
multiplier and  Tektronix  Sampling  Oscillosco,  No  pulse  broadening  was 
evident  and  it  was  assumed  that  the  magnitude  of  the  broadening  was  beyond 
the  resolution  capability  of  this  measurement  system. 

The  measurement  technique  adopted  is  analogous  to  that  used  to  measure 
group  delay  in  electrical  transmission  lines.  The  application  of  this 
method  requires  a periodic  radio  frequency  (RF)  signal,  a condition  satis- 
fied by  the  mode  locked  lie:Ne  laser.  The  incident  and  scattered  intensities 
are  used  to  derive  two  RF  signals  whose  relative  phase  re.' ationships 
are  compared.  Time  delay  is  then  deduced  from  the  phase  measurements. 

C.  Analysis  of  measurement  errors 

The  only  measurement  errors  considered  are  those  relating  to  the 
relative  phase  between  incident  and  scattered  signal  waveforms.  These 
errors  may  result  from 
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1.  Difference  in  relative  signal  levels. 


2.  Systematic  variation,  or 

3.  Detector  translation 

1.  Difference  in  relative  signal  levels 

The  primary  range  of  this  error  occurs  when  Channel  B amplitude 
varies  between  5 mV  and  25  mV  producing  a phase  change  of  about  2.5  degrees. 
Below  1 mV  the  phase  changes  due  to  changing  relative  amplitude  are 
negligible. 

A wedge  attenuator  with  10:1  intensity  attenuation  ratio  was  used 
to  obtain  a calibration  curve  (phase  correction  vs  Channel  B amplitude) 
and  all  raw  data  was  adjusted  using  the  appropriate  correction  prior  to 
analysis. 

2.  Systematic  variation 

In  general,  the  phase  fluctuated  about  a mean  value  during  measure- 
ment with  average  variation  of  +0.2  degrees.  The  phase  values  recorded 
were  obtained  by  estimating  the  mean  value  during  the  observation  interval 
usually  5 to  10  seconds.  Additionally,  a chart  recorder  was  used  on 
selected  data  runs  to  illustrate  typical  variation  conditions  (Fig.  2-4). 

3.  Detector  translation 


Mechanical  slack  at  the  swing  arm  hinge  point  allows  slight  transla- 
tion of  the  detector.  However,  the  maximum  variation  is  0.01mm  and  the 


resulting  phase  error  is  negligible. 


D.  Summary  of  results 


1.  Effect  of  experiment  parameters  on  phase  shift 


a.  Detector  polarizer 

Measurements  were  made  with  the  polarizer  oriented  both  parallel  and 
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(arbitrary  units)  Aa  (arbitrary 


perpendicular  to  the  incident  beam  polarization.  No  significant  difference 
in  the  phase  shift  was  apparent. 

b.  Aperture 

Several  apperture  sizes  were  used.  The  phase  shift  and  aperture  size 
are  directly  related  since  the  aperture  controls  the  FOV. 

c.  Scatterer  density 

Several  different  concentrations  were  used,  each  producing  similiar 

phase  shifts.  At  the  highest  density  the  signal  level  at  the  vector 

voltmeter  was  insuffi  lent  for  reliable  and  repeatable  measurements,  thus, 
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the  scatterer  density  wa  restricted  to  a rather  narrow  range  between  10  cm 
and  A x lO^cm 

2.  Review  of  results 

The  results  were  obtained  with  detector  aperture  and  FCV  as  described 
in  Section  l.C.  The  minimum  phase  shift  measured'  was  0.7  degrees;  the 
maximum  was  1.9  degrees  and  the  mean  was  1.25  degrees  for  21  trials.  At 
low  scatterer  densities  the  intensity  of  the  central  beam  caused  excessive 
photomultiplier  current,  therefore,  all  phase  measurements  were  taken 
between  <j>  = 2°  and  <|>  = 12°  (4>  = 2°  corresponds  to  a detector  position 
clear  of  the  central  beam).  A typical  plot  is  shown  in  Fig.  2-5. 

At  the  maximum  scatterer  density  the  effective  optical  thickness  of 
the  medium  was  slightly  less  than  2.  Thus,  i'c  is  not  surprising  that  the 
phase  shift  is  comparable  to  that  expected  from  single  scattering  since  the 
expectation  of  the  number  of  scattering  events  for  a given  photon  is  on 
the  order  of  1 or  2. 


CONCLUSIONS 


The  tine,  delays  measured  in  this  experiment  did  not  produce  signifi- 
cant time  spreading  of  the  He:Ne  mode  locked  pulse.  Neither  would  they 
produce  significant  time  spreading  of  much  shorter  pulses  typical  of 
operational  optical  coomunication  systems.  This  is  consistent  with  the 
fact  that  single  scattering  was  the  dominant  mechanism  producing  the 
nultipath  delays.  It  has  been  pointed  out  that  the  relatively  low 
intensity  (25  mW)  of  the  optical  source  restricted  the  maximum  scattering 
medium  density  which  could  be  used  while  assuring  adequate  intensity  at 
the  detector.  This  restriction  lead  to  predominantly  single  scattering. 

The  measurement  of  multipath  delay  by  using  the  comparison  of  relative 
phase  of  the  (periodic)  incident  and  scattered  pulse  envelopes  circumvents 
the  detector  rise  time  limitations  and  could  be  applied  tc  measurements 
with  much  shorter  pulses. 

Further  experimentation  is  needed  using  shorter,  more  intense 
pulses.  One  possibility  would  be  a mode  locked  Argon  ion  laser.  Such 
a source  could  provide  pulses  in  the  300  picosecond  range  and  power  on 
the  order  of  watts.  The  increased  intensity  is  essential  to  increasing 
the  optical  thickness  of  the  medium  producing  multiple  scattering  evenls. 

For  large  optical  thickness,  Bucher’s  computer  simulations  predict 
multipath  delays  of  nanoseconds.  If  the  main  beam  has  suffered  great 
extinction  and  the  receiver  field  of  view  has  been  increased  to  collect 
scattered  photons,  the  delayed  photons  may  produce  very  significant  time 
spreading  of  the  received  pulse.  For  example,  a 300  picosecond  pulse 
system  with  300  picosecond  guard  time  has  a potential  data  rate  of 
1.67  Gigabits/sec.  £n  the  event  the  pulses  are  spread  by  multipath  pro- 
cesses to  1 nanosecond  the  data  rate  is  reduced  to  770  Megahits/sec. 
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The  effects  of  multipath  pulse  spreading  on  the  potential  data  rate 
of  currently  feasible,  state  of  the  art  optical  communications  systems 
are  far  more  profound.  Consider  a mode  locked  dye  laser  capable  >f 
generating  1 picosecond  pulses;  beam  superposition  techniques  can  be 
used  to  produce  a periodic  pulsed  output  with  period  of  several  picoseconds 
(say  10  picoseconds);  in  a synchronized  link,  such  pulses  can  be  resolved 
with  optoelectronic  sampling  [10].  Such  a system  has  a potential  data 
rate  of  100  Gigabits/sec;  however,  multipath  pulse  spreading  of  1 nano- 
second would  reduce  it  by  a factor  of  100.  Thus  it  is  apparent  that 
multipath  pulse  spreading  produces  severe  limits  for  ultra  short  pulse, 
high  data  rate  optical  communications  systems. 
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APPENDIX  A 

Let  Ij_  be  defined  as  the  Intensity  of  the  wave  scattered  by  a single 
sphere  from  an  incident  beam  of  unit  intensity  with  polarization  per- 
pendicular to  the  observation  plane. 


! -_iV  £l 

A 4n2r2  ? 


(A.  1) 


where  A is  the  wavelength  in  the  medium  and 


a3  = *3  . (iM)3 


(A.  2) 


as  previously  defined  and  — r is  the  tabulated  value  for  the  far  field 

a 


(r  » A)  at  the  angle  y. 

Since  y is  the  angle  of  observation  with  respect  to  the  reverse 
direction  of  the  incident  beam,  y is  related  to  9 by 


0 = 7i  - y 


(A. 3) 


The  following  values  for  ~ as  a function  of  Y assume  a = 3.2,  m - 1.2, 

or 

(Table  A-l): 
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O = 3.2 


m=  L20 


r 

b./a3 

/ 

in  /a3 

h/a* 

0 

0.0014914 

0.0014914 

0.0029828 

5 

0.0014827 

0.0014641 

0.0029468 

10 

0.0014655 

0.0013861 

0.0028516 

15 

0.0014650 

0.0012687 

0.0027337 

20 

0.0015198 

0.0011277 

0.0026475 

25 

0.0016770 

0.00097987 

0.0026569 

So 

0.0019853 

0.00083922 

0.0028246 

35 

0.0024863 

0.00071441 

0.0032007 

40 

0.0032041 

0.00060789 

0.0038119 

45 

0.0041388 

0.00051714 

0.0046559 

50 

0.0052524 

0.00043922 

0.0056916 

55 

0.0064624 

0.00037582 

0.0068382 

60 

0.0076453 

0.C0033728 

0.0079826 

65 

0.0086325 

0.00034508 

0.0089775 

70 

0.0092297 

0.00042743 

• 0.0096571 

75 

0.0092484 

0.00061006 

0.0098584 

80 

0.0085504 

0.00090357 

0.0094540 

85 

0.0071135 

0.0012955 

0.0084090 

90 

0.0050997 

0.0017559 

0.0068556 

95 

0.0029466 

0.0022751 

0.0052217 

100 

0.0014406 

0.0029429 

0.0043836 

105 

0.0017651 

0.0040731 

0.0058382 

110 

0.0055021 

0.0063617 

0.011864 

115 

0.014566 

0.011042 

0.025608 

120 

0.031043 

0.019976 

0.051019 

125 

0.056955 

0.035614 

0.092569 

130 

0.093922 

0.060756 

0.15468 

135 

0.14281 

0.098078 

0.24089 

140 

0.20339 

0.14949 

0.35288 

145 

0.27406 

0.21539 

0.48944 

150 

0.35173 

0.29398 

0.64570 

155 

0.43195 

0.38097 

0.81292 

160 

0.50925 

0.46973 

0.97899 

165 

0.57767 

0.55193 

1.1296 

170 

0.63145 

0.61873 

1.2502 

175 

0.66585 

0.66248 

1.3283 

180 

0.67769 

0.67769 

1.3554 

Tabic  A-l 


s^. 


■ifiZ'MzMt  ; 


39 

APPENDIX  B 

Determination  of  N 
o 

The  number  of  particles  per  unit  volume,  N was  determined  directly 
by  calculation  from  specifications  provided  by  the  manufacturer  and 
indirectly  by  measurement  of  the  extinciton  coefficient-  The  values  ob- 
tained in  either  case  agree  closely. 

Calculation  of  N 

The  weight  of  a single  latex  sphere,  W = SV  where  5 is  the  density 
(of  polystyrene  = 1.05  g/ml)  and  V is  the  volume  of  the  sphere.  Since 
V = where  D is  the  diameter  = 0.481  pm,  then  W = (1.05  g/ml)^n 

__  C *1  — 1 / 

x (2.405  x 10  cm)  = 6.12  x 10  g.  Assume  a 1 ml  sample  of  the  latex 
particles  in  aqueous  suspension  weighs  1 g.  The  manufacturer  asserts  that 
the  sample  contains  10%  solids  by  weight.  Then  the  number  of  particles 
in  a 1 ml  sample,  N is  the  total  weight  of  particles  divided  by  the 
weight  of  one  particle. 


N = 


= 1.63  x io12 
6.12  x io"14  g 


(B.l) 


Measurement  of  N 

The  extinction  coefficient,  £ was  determined  by  measuring  the  trans- 
mitted intensity  of  the  He:Ne  laser  main  beam  after  propagating  through 

a cell  containing  a highly  dilute  suspension  of  the  latex  particles  in 

2 

deionized  water.  Using  the  relation  I = l^exp(-^L),  £ = Nna  Qgxt 
where 

- transmitted  vntvr:-ity 
1^  = incident  intensity 
I.  * cell  length 
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number  of  particles 
radius  of  particle 


extinction  efficiency  =0.7492 
(given  by  Wickramasinghe  [HJ) 


Solving  for  N we  find 


N = 1.94  x 10 


(B.2) 


Calculation  of  N 


The  concentration  of  scatterers  in  the  chamber  per  ml  of  latex 


particles  sample,  is  equal  to  — where  V is  the  chamber  volume  = 
3000  ml.  Using  the  measured  value  for  N 


„ 1.94  x 10**  ln8  -3 

N = r = 6.46  x 10  cm 

3 x 10  ml 


(B.3) 


A plot  of  Nq  versus  volume  of  latex  particles  is  given-  in  Fig.  B-l. 
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APPENDIX  C 

Calculation  of  Optical  Thickness 
The  optical  thickness,  T is  defined  as 


It  I 

L I 


T - T/D 


(C.l) 


where  T is  the  physical  thickness  of  the  medium  and  D is  the  mean-free 
path  between  scatters.  Equation  (C.l)  is  valid  for  isotropic  scatter; 
for  the  case  of  nonisotropic  scatter  the  mean-free  path  may  be  replaced 
by  the  diffusion  distance  [12J 


®d  1 - <cos  9> 


(C.2) 


si 

Is 


i . 


where  <cos  0>  is  the  average  cosine  of  the  scattering  angle  relative  to 
the  incident  direction. 

In  a homogeneous  body  of  scatterers,  D is  the  Inverse  of  the 
extinction  coefficient,  £.  If  the  scatterers'  radii  are  uniform 


D = — « 

5 Q N ira2 
ext  o 


(C.3) 


where  Q Is  the  extinction  efficiency  factor,  NQ  is  the  number  of 
scatterers  per  unit  volume,  and  a is  the  scatterer's  radius.  Substituting 
(C.3)  into  (C.2)  we  obtain: 


2 M 

Q jra  (1  - <cos  0>)  o 
ext  ' 


(C.4) 


Prom  (C.l)  and  (C.2)  wo  obtain  an  effective  optical  thickness 


T Q . ira  (1  - <cos  0>)-N 
ext  o 


(C.5) 


Solving  (C.4)  and  (C.5)  for  and  (given  that  T = 25.4  cm,  <cos  0> 
* 0.9,  Qext  “ 0.7492  [10],  and  a = 0.481  pm)  we  find 


Dd  = 7.345  x 109  cm**2  * ^ 

o 

T.  - 3.458  x l(f9  cm3  * N 
d o 


(C.6) 

(C.7) 


Plots  of  D and  T versus  N are  given  in  Fig.  C-l. 


Fig.  C-l.  Diffusion  distance  and  optical  thickness 
T.  versus  concentration  N . 


